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We study the nature of the energy transfer process within a pair of coupled two-level systems 
(donor and acceptor) subject to interactions with the surrounding environment. Going beyond a 
standard weak-coupling approach, we derive a master equation within the polaron representation 
that allows for investigation of both weak and strong system-bath couplings, as well as reliable 
interpolation between these two limits. With this theory, we are then able to explore both coherent 
and incoherent regimes of energy transfer within the donor-acceptor pair. We elucidate how the 
degree of correlation in the donor and acceptor fluctuations, the donor-acceptor energy mismatch, 
and the range of the environment frequency distribution impact upon the energy transfer dynamics. 
In the resonant case (no energy mismatch) we describe in detail how a crossover from coherent 
to incoherent transfer dynamics occurs with increasing temperature [A. Nazir, Phys. Rev. Lett. 
103, 146404 (2009)], and we also explore how fluctuation correlations are able to protect coher- 
ence in the energy transfer process. We show that a strict crossover criterion is harder to define 
when off-resonance, though we find qualitatively similar population dynamics to the resonant case 
with increasing temperature, while the amplitude of coherent population oscillations also becomes 
suppressed with growing site energy mismatch. 
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I. INTRODUCTION 

A fascinating series of recent experiments demonstrat- 
ing signatures of quantum coherence in the energy trans- 
fer dynamics of a variety of systems 1 -"— has sparked re- 
newed interest in modeling excitation energy transfer 
beyond standard methods^"— This process, which oc- 
curs when energy absorbed at one site (the donor) is 
transferred to another nearby site (the acceptor) via a 
virtual photon,— is often considered to be incoherent; 
the result of weak donor-acceptor interactions, treated 
perturbatively using Fermi's golden rulei 21 ' 22 However, 
though this approach has proved to be immensely suc- 
cessful when applied in many situations , 23 ! 24 accounting 
for quantum coherence within the energy transfer dynam- 
ics requires an analysis beyond straightforward perturba- 
tion theory in the donor-acceptor interaction. 

An alternative starting point for investigations into co- 
herent energy transfer is to treat the system-environment 
interaction as a perturbation instead. Such weak- 
coupling theories, often referred to as being of Redfield 
or Lindblad type depending upon the approximations 
made in their derivation,— have been successfully ap- 
plied to elucidate a number of effects that could be at 
play in multi-site donor-acceptor complexes. Examples 
include studying the interplay of coherent dynamics and 
dephasing in promoting efficient energy transfer in quan- 
tum aggregates*^"— exploring the role of environmen- 
tal correlations in tuning the energy transfer process,- 3 " 3 - 
and extensions to assess the potential importance of non- 
Markovian effects^^ 



Nevertheless, in order to properly understand the tran- 
sition from coherent to incoherent energy transfer which 
occurs as the system-environment coupling or tempera- 
ture is increased ) 11 ! 36 " — it is necessary to be able to de- 
scribe the system dynamics beyond either of these limit- 
ing cases . 24 ! 39 ' 40 Building on earlier work , 36 ' 41 " — a num- 
ber of methods have been put forward to accomplish 
this. For example, modifications to both Redfield 4 ^"— 
and Forster— "— theory have extended the range of valid- 
ity of both approaches. Moreover, it is possible to define a 
new perturbation term through the small polaron trans- 
formation^ which under certain conditions allows in- 
terpolation between the Redfield and Forster limits.—"— 
For particular forms of system-environment interaction, 
this can also be achieved through the hierarchical equa- 
tions of motion techniquej 16 ' 52 Numerically exact calcula- 
tions, based, for example, on path integral , 15 ' 53 numerical 
renormalisation group 5 ^ and density matrix renormalisa- 
tion groupi 4 - methods, have also been applied to study 
energy transfer beyond perturbative approaches. 

In this work, we investigate the conditions under which 
coherent or incoherent motion is expected to dominate 
the energy transfer dynamics of a model donor-acceptor 
pair. Following Ref. [TJ we employ a Markovian mas- 
ter equation derived within the polaron representation 
for this purpose, since it allows for a consistent analysis 
of the dynamics from weak to strong system-bath cou- 
pling (or, equivalently, low to high temperatures) j 5 1 ' 55 
In addition to presenting a full derivation of the theory, 
we also extend it to explore in detail the important ef- 
fects of donor- acceptor energy mismatch, deriving ana- 
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lytical forms for the dissipative dynamics valid over a 
large range of parameter space. Furthermore, we move 
beyond the scaling limit studied in Ref. [ll| to consider an 
environment frequency distribution of finite extent, char- 
acterised by a high-frequency cut-off in the bath spectral 
density. In the resonant case (no energy mismatch) we 
define a strict crossover temperature above which the en- 
ergy transfer dynamics ceases to be coherent 

Of particular practical interest is the role played by 
correlations between the donor and acceptor environmen- 
tal fluctuations, suggested as one mechanism by which 
quantum coherence may survive in the energy transfer 
process under otherwise adverse conditions^ & 11 i 53 ' 56 ~ — 
(though see, for example, Refs. |6l| and [62| for alterna- 
tives). These correlations are also easily treated within 
our formalism, through position-dependent couplings be- 
tween the system and the common environment. As the 
donor and acceptor are brought closer together, there 
comes a point at which their separation becomes com- 
parable to, or smaller than, the wavelength of relevant 
modes in the bath. As this happens, fluctuations at each 
site become ever more correlated, and dephasing effects 
are suppressed. We shall show, consequently, that as the 
level of correlation increases, so too does the crossover 
temperature to the incoherent regime. Hence, strong cor- 
relations lead to the survival of coherence at high tem- 
peratures. 

Off-resonance, we find that it is less straightforward to 
define a crossover temperature. In contrast to the res- 
onant case, for sufficient energy mismatch between the 
donor and acceptor, increasing the temperature causes 
the amplitude of the coherent contribution to decrease, 
though not to disappear altogether. In principle, it then 
becomes possible for a coherent component to exist in 
the dynamics at all but infinite temperatures. Although 
we are then unable to define a crossover in quite the same 
way, we still find that bath correlations have a qualita- 
tively similar effect to the resonant case, protecting co- 
herence in the transfer process. 

The paper is organised as follows. In Section [TT| we 
introduce our model, and derive a master equation de- 
scribing the donor-acceptor dynamics within the polaron 
representation. Section IIIII considers the resonant case 
and the coherent-incoherent crossover. In Section [IVI we 
investigate off-resonant energy transfer and obtain ana- 
lytic expressions for the dynamics in a number of limits. 
Finally, in Section [V] we summarise our results. 



II. POLARON TRANSFORM MASTER 
EQUATION 

A. The system and polaron transformation 

We consider a donor-acceptor pair (j = 1,2), each site 
of which is modeled as a two-level system with ground 
state \G)j, excited state \X)j, and energy splitting Cj, 
The pair interact via Coulombic energy transfer with 



strength V, which is responsible for the transfer of exci- 
tation from one site to the other. We label the state cor- 
responding to a single excitation on site 1 as |1) = \XG), 
and that on site 2 as |2) = \GX). The environment sur- 
rounding the donor-acceptor pair is modelled as a com- 
mon bath of harmonic oscillators, coupled linearly to the 
excited state of each site. The total system-bath Hamil- 
tonian in the single excitation subspace (in which energy 
transfer occurs) is therefore written (where h = 1) 

ffsuB = e 1 |l)<l|+e 2 |2X2| + n|l><2| + |2><l|) 
+ |l)(l|i?( 1 ) + |2)(2| J B( 2 )+^ Wk 6j i 6 k , 

k 

(1) 

where the bath is described by creation (annihilation) 
operators (6 k ) with corresponding angular frequency 
Wk, and wavevector k. The bath operators are given by 

5z = X^klSk ^k+S'k *&k)i w ith coupling constants g^K 
As in Ref.llll we shall consider the case in which each site 
is coupled to the bosonic bath with the same magnitude 
| (7k |) but make the separation between the sites explicit 
through position-dependent phases in the coupling con- 
stants of the form = \gk\e lk ' Tj , with Vj being the 
position of site j. As we shall see, this form of coupling 
gives rise to correlations between the bath influences ex- 
perienced at each site, allowing a range of totally corre- 
lated, partially correlated, and completely uncorrelated 
fluctuations to be explored] 11 i 53 

A standard weak-coupling approach to the system dy- 
namics would now be to derive a master equation for the 
evolution of the reduced system density operator under 
the assumption that the system-bath interaction terms, 
as written in Eq. (fTJ), can be treated as weak perturba- 
tionsj 25 i 63 In this work, we shall instead derive a master 
equation describing the donor-acceptor energy transfer 
dynamics in the (now widely used) polaron representa- 
tion^i ~ 13 i 43 ' 51 i 55 ' 64 whereby we displace the bath oscil- 
lators depending on the system state. We may then 
identify alternative perturbation terms, which can be 
small over a much larger range of parameter space than 
those in the original representation. In particular, the 
polaron framework allows us to reliably explore from 
weak (single-phonon) to strong (multiphonon) coupling 
regimes between the system and the bath, provided that 
the energy transfer interaction V does not become the 
largest energy scale in the problem (in which case the 
full polaron displacement is no longer appropriate^) , 
and that there is no infra-red divergence in its bath- 
renormalised value Vr (see Eq. J7| below) ^ In contrast, 
with a weak system-bath coupling treatment we would 
only be able to probe single-phonon bath-induced pro- 
cesses, and hence not be able to properly explore the 
crossover from coherent to incoherent dynamics in which 
we are primarily interested. 

To proceed, we thus apply a unitary transformation 
which displaces the bath oscillators according to the lo- 
cation of the excitation. Defining Hp — e s i7suBe~ s , 
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where 



S = |l)<l|PG£Vwk) + \2)(2\ P(g^/u; k ), (2) 



with bath operators P(a^) — J^kC ^!— o^k)) results in 
the polaron transformed spin-boson Hamiltoniar* 5 -^ Hp = 
H +Hi, with 



and 



Hi = V(B x a x 



ByCXy) 



(3) 



(4) 



Here, the bias e = t\ — €2, gives the energy difference 
between the donor and acceptor, while the Pauli opera- 
tors are defined in a basis in which a z = |1)(1| — \2)(2\ = 
\XG)(XG\ - \GX){GX\. The bath operators appearing 
in Eq. (0 are constructed as B x = (1/2)(B + + B- - 2B) 
and B v = (i/2)(B+ — B-), where 



b ± = Y[d(± 




(5) 



with displacement operators D(±au) = exp[±(ak& k — 
a k 6k)]- Note that the interaction terms in Eq. J2J there- 
fore depend upon the difference in donor and acceptor 
system- bath couplings an d respectively. Im- 

portantly, the term driving coherent energy transfer in 
Eq. © will not be treated perturbatively, though it does 
now have a bath-renormalised strength, Vr = BV, where 



B 



exp 



E 

k 



9k 



(1 - cos(k ■ d)) coth(/3w k /2) 



(6) 



is the expectation value of the bath operators with re- 
spect to the free Hamiltonian: B = (B±)h - The donor- 
acceptor separation is given by d = ri — r 2 . 

In order to calculate the renormalisation factor, we 
take the continuum limit to convert the summation in 
Eq. ([5]) into an integral. Defining the bath spectral den- 
sity J(oS) — Ylk l3k| 2 <5(^ ^ w k ), which contains infor- 
mation regarding both the density of oscillators in the 
bath with a given frequency, and also how strongly those 
oscillators interact with the donor-acceptor pair, and as- 
suming a linear, isotropic dispersion relation, we find 



B 



exp 



£M (1 



F D (u,d))coth(j3ui/2) 



(7) 



Here, /3 = 1/fcgT is the inverse temperature, while the 
function Fu(uj,d) captures the degree of spatial corre- 
lation in the bath fluctuations seen at each site, and 
is dependent upon the dimensionality of the system- 
bath interaction (D = 1,2,3).—^^ In one dimen- 
sion Fi(u>,d) = cos(ud/c), with c the bosonic exci- 
tation speed, in two dimensions F2(u>,d) — Jo(u>d/c), 
where Jq(x) is a Bessel function of the first kind, and 



in three dimensions F%(w,d) = sinc(u)d/c). In all cases 
Fd(w, d) — > I as d — > 0, i.e. when the donor and accep- 
tor are at the same position, bath fluctuations are per- 
fectly correlated, and the energy transfer strength is not 
renormalised (Vr — > V). In fact, in this limit dissipativc 

process are entirely suppressed (provided \g k \ — \) 
and energy transfer remains coherent for all times and in 
all parameter regimes in our model (the single-excitation 
subspace is then decoherence-free i 67 i 68 ) In two and three 
dimensions, as d — > 00, Fjj(uj,d) —¥ 0, and the renor- 
malisation takes on the value that would be obtained 
by considering separate, completely uncorrelated baths 
surrounding the donor and acceptor. In the following, 
we shall characterise the degree of correlation in terms 
of the dimensionless parameter fi = c/u>od, where luq is 
a typical bath frequency scale (see Eq. (|2T))) below). We 
therefore have fi = in the absence of correlations, /1 < I 
for weak correlations, and \x > 1 for strong correlations. 



B. Markovian master equation 

Having identified a new perturbation term by trans- 
forming our Hamiltonian to the polaron representation, 
we can now construct a master equation describing the 
evolution of the donor-acceptor pair reduced density op- 
erator p up to second order in Hj. We employ a standard 
Born-Markov approach, which yields a polaron frame, in- 
teraction picture master equation of the form 25 



d~ P (t) 
dt 



drttB [[Hi(t) 7 [Hi(t-T),p(t)®p B ]] , (8) 



where tildes indicate operators in the interaction picture, 
0(t) = e lH ° t Oe~' lHot , and trs denotes a trace over the 
bath degrees of freedom. In deriving Eq. ((SJ we have 
assumed: (i) factorising initial conditions for the joint 
system-bath density operator within the polaron frame, 
X(0) = p(0) ® p B , with p B = e-' mB /tr B (e~ imB ) be- 
ing a thermal equilibrium state of the bath; (ii) that 
by construction the interaction is weak in the polaron 
frame so that we may factorisc the joint density operator 
as x(t) = p(t)pB at all times; (iii) that the timescale 
on which the donor-acceptor system evolves apprecia- 
bly in both the Schrodinger and interaction pictures is 
large compared to the bath memory time tb- Since, for 
the spectral density we shall consider below, tb ~ l/w c , 
where u> c is a high-frequency cutoff (see Eq. ([2"5|) ). this 
is not too restrictive, as we must keep V < uj c any- 
way in order for the polaron theory to work well. We 
note that interesting non-Markovian and non-equilibrium 
bath effects have been explored in the polaron formalism 
in Refs.lllandlll 

Inserting Eq. Q into Eq. ([8]) , and moving back into the 
Schrodinger picture, we arrive at our Markovian master 
equation describing the energy transfer dynamics within 
the single-excitation subspace, and written in the polaron 
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frame as 
dp® 



Of 



= -i[(e/2)a z + V R a x ,p(t)} 



-V 2 dr([CT x , a x (-r)p(t)]A xx (r) 

+ [(Jy, <Jy\ — T)p(t)}A yy (r)+R.c.), 

(9) 

where H.c. denotes Hermitian conjugation. The effect 
of the bath is now contained within the correlation func- 
tions Au(t) — (Bi(t)Bi(0))h o , which are given explicitly 
by 



A xx (t) = (5 2 /2)(e 0(T) + e^M - 2), 
A yy (r) = (B 2 /2)(e <AW -e-^)), 



where 



f c 

(r) = 2 / 
Jo 



did 



^(l-F D (u,d)) 



x (coswr coth(f3u>/2) — i sinwr) 



(10) 
(11) 



(12) 



Notice that the phonon propagator, </i(r), is correlation- 
dependent due to the factor (1 — Fp,(ui, d)), and so clearly 
the dissipative effect of the bath will be dependent upon 
the degree of correlation too. For example, as d — > 0, 
Fd(u>, d) — > 1, and the dissipative contribution to Eq. © 
vanishes, as anticipated earlier. 



C. Evolution of the Bloch vector 

We solve our master equation in terms of the Bloch vec- 
tor, defined as a — (a x ,a y ,a z ) T = ((a x ) , (a y ) , (a z )) T . 
As we are working exclusively in the single-excitation 
subspace, a x and a y describe the coherences between the 
states |1) = \XG) and |2) = \GX), while a z captures the 
donor-acceptor population transfer dynamics generated 
by the coupling V. 

Though Eq. Q is written in the Schodinger picture, 
it is still in the polaron frame, and so we must de- 
termine how expectation values in the polaron frame 
are related to those in the original, or "lab" frame. 
We can see this by writing a,; = tr s+b (fiXi(i)) = 
tr s+ B{cne~ s x(t)e s ) = tr s+B {e s aie~ s p{t)p B ), where 
Xl(0 — e ~ S x{t) eS is the lab frame total density op- 
erator, and we have made use of the Born approxima- 
tion in the polaron frame to write x(t) = p{t)PB- Since 
e s a x e- s = |2)(1| S_ + |l)(2| B+, e s a y e~ s = i(\2)(l\ B_- 
|1)(2| B+), and e s a z e~ s — a z , this implies that the lab 
Bloch vector elements are a% — Baip, for i = x,y, and 
ol z = a z p, where oup is an expectation value in the po- 
laron frame: c^p = Tr,g((7jp(t)). Alternatively, we can 
define a matrix L which maps the polaron frame Bloch 
vector {dp) to its lab frame counterpart (a): a = L-ap, 
where L — diag(_B, £?, 1). 



Working in terms of the Bloch vector, we arrive at an 
equation of motion of the form 



a(t) = M ■ a(t) + b. 



(13) 



In the following, we shall often be interested in deter- 
mining whether the energy transfer dynamics is predom- 
inantly coherent or incoherent. It is then helpful to write 
Eq. as 



a'(t) = M ■ a'(t), 



(14) 



with tt'(i) = ct(t) — ck(oo), where a(oo) = — M^ 1 ■ b is 
the steady state. This makes clear that the nature of the 
energy transfer process lies solely in the matrix M, while 
the inhomogeneous term b is needed only in determining 
the steady state. 

Equipped with the eigensystem of M, we may deter- 
mine the corresponding time evolution as follows: an 
eigenvector of M, say m^, has equation of motion rn L — 
qiirii, where qi is the corresponding eigenvalue. Its sub- 
sequent evolution then has the simple exponential form 
rrii{t) — rriie qit . More generally, we can say that any 
initial state ct'(0) will have subsequent evolution 



a'(t) — aiiriie' 



<i,i 



(15) 



where the coefficients ai are determined by the initial 
conditions (i.e. the solutions of a'(0) = J2i a i rn i)- The 
solution to the full inhomogeneous equation is then found 
simply by addition of the steady state: ait) — a'(t) + 
a(oo). 



III. RESONANT ENERGY TRANSFER 

We start by considering the important special case of 
resonant energy transfer, in which the interplay of co- 
herent and incoherent effects is particularly pronounced. 
As we shall see, in this situation it is relatively straight- 
forward to derive a strict criterion governing whether or 
not we expect the energy transfer dynamics to be able to 
display signatures of coherencei 11 ' 51 Hence, resonant con- 
ditions provide a natural situation in which to begin to 
understand, for example, the role of bath spatial correla- 
tions^ i 11 ' 33 ' 53 ' 56 " — or the range of the bath frequency dis- 
tribution in determining the nature of the energy transfer 
process. 

Setting the donor-acceptor energy mismatch to zero, 
e = 0, wc find from Eq. @ dynamics generated by an 
expression of the form a — Mr ■ a + bp , with 



M R = 





B~ 1 {2V R 





-2BV R 

As) -r 4 



(16) 
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and b R = (— Bk x , 0, 0) T , where 
T y = 2V 2 lxx (0), 



(17) 

= V 2 ( lyy {2V R ) + lyy (-2V R )) + 2V 2 lxx (0), (18) 

A 3 = 2V 2 (Syy(2V R ) - Syy(-2V R )), (19) 

«x = V*(7 w (2Vr) - J yy (-2V R )). (20) 

The rates and energy shifts are related to the response 
functions 

K H {u))=J^ &Te WT k H {T) = ^ llt {u:)+iS lt {u:), (21) 



such that 

7«(w) = 2Re[JT«(w)] 



+ 00 



dre~ r A«(r), (22) 



and Su{w) = Im[ii'fi(a;)]. 

The resonant steady-state is straightforwardly found 
to be 



a x (co) = -Btanh(0V R ), 



(23) 



while a y (oo) = a z (oo) = 0. Notice that while this is sim- 
ilar in form to the steady-state that would be obtained 
from a weak system-bath coupling treatment) 25 ! 69 ce x (oo) 
is determined here by V R , rather than the original cou- 
pling V, and there is also an extra factor of B suppressing 
its magnitude. 

The eigenvalues of M R are given by qi = T y — T z and 
92 = q* 3 = -(l/2)(r„ + T Z + Thus, referring to 

Eq. (fl5|) . we see that 



^8V r (2V r + \ 3 ) - (r 2 -T y ) 2 



(24) 



determines whether or not any coherence exists within 
the energy transfer dynamics. Considering the initial 
state a(0) = (0, 0, 1) T , corresponding to excitation of 
the donor, p(0) = |1)(1| = \XG}(XG\, we find analytical 
forms for the evolution of the Bloch vector components: 



a x {t) 

a y {t) = - 
a z (t) 



-Bt&nh{/3V R ){1 
2BV R 



e - ( r,+r,) t /2 sin 



= c -(r H +r,)t/2 



cos 



(¥) 



2 

r„ 



(25) 
(26) 



sin 



(¥) 



(27) 



Inspection of Eqs. (f2~4"|) and (|2T|) allows us to identify a 
crossover from coherent to incoherent motion in the en- 
ergy transfer dynamics as the point at which oscillations 
in the population difference vanishi^i 



(T, 



8V R (2V R + X 3 ). 



(28) 



For (T z - Ty) 2 < 8V R (2V R + A 3 ), £ R is real and both the 
population difference and coherence a y describe damped 
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FIG. 1: Population difference as a function of scaled time 
wot for temperatures of ksT /ujq — 1 (blue dashed curve), 
UbT /uo = 5 (green dotted curve), ksT/uJo = 12 (orange solid 
curve) and fcsT/cjo = 20 (red dot-dashed curve). Parameters: 
a = 0.05, V/0J0 = 0.5, uj c /ljq — 4, e = and /i = c/ljqcI = 0.5. 



oscillations, while for (T z - Y y ) 2 > 8V R (2V R + A 3 ), £ R 
is either zero or imaginary, with the resulting dynamics 
then being entirely incoherent. 

To further analyse the behaviour of a z (t), and the con- 
ditions for which the boundary defined by Eq. (|28p is 
crossed, we now take a specific form for the system-bath 
spectral density. For a large enough bath we may approx- 
imate J(uj) = ^ k |.9k| 2 ^(^k — w) as a smooth function of 
uj. In this work we consider a spectral density of the form 



J(w) = a— T-e 



u/a 



(29) 



where a is a dimensionless quantity capturing the 
strength of the system-bath interaction, and ujq is a typi- 
cal frequency of bosons in the bath, which sets an overall 
energy scale. The cubic frequency dependence in Eq. (|29|) 
is typical, for example, in describing dephasing due to 
coupling to acoustic phononS ) 51 ' 70 but can also be used 
to elucidate the behaviour in which we are interested in 
general^ The cut-off frequency uj c is needed to ensure 
that vacuum contributions remain finite, and is related 
to parameters specific to the particular physical system 
one wishes to model. The inverse cut-off frequency also 
sets a typical relaxation timescale for the bathi^ 5 . 

To illustrate the dynamics and crossover behaviour in 
the resonant case, in Fig. [1] we plot the population differ- 
ence (a z ) as a function of the scaled time oj^t for a range 
of temperatures, showing the transition from coherent 
to incoherent transfer as the temperature is increased. 
In this plot, and all the following, we consider three- 
dimensional coupling, F 3 (uj,d) — sinc(cjG?/c). The role of 
bath spatial correlations in protecting coherence can be 
seen in Fig. [21 where we again plot the evolution of the 
population difference (the insets show the correspond- 
ing coherence a y ), this time for representative interme- 
diate and high-temperature cases. The different plots in 
Fig. [21 correspond to zero correlations, characterised by 
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FIG. 2: Population difference as a function of scaled time 
coot for temperatures of feT/uo = 5 (blue dashed curves) 
and ksT/ujo = 10 (red dotted curves), and for separations 
corresponding to no correlation, \i = c/u>od = (top), weak 
correlations, /j = 0.5 (middle), and strong correlations fi = 2 
(bottom). The insets show the evolution of the correspond- 
ing coherence a y . Parameters: a = 0.05, V/cjo = 0.5, and 
lo c /loo — 4. 



/i = c/u>od = (d — > oo, top), weak correlations, /j = 0.5 
(middle), and strong correlations /i = 2 (botttom).— Pro- 
gressing from the uppermost plot to the lowest, we clearly 
see that an increase in correlation strength prolongs the 
timescale over which oscillations in both the population 
difference and coherence persist. Moreover, by looking at 
the curves corresponding to the higher temperature (red, 
dotted), we can see that as the degree of correlation is 
increased from zero, the dynamics moves from a regime 
showing purely incoherent relaxation, to a regime which 
displays coherent oscillations at the same temperature. 



The increase in correlations is thus able to extend the 
region of parameter space which permits coherence^ as 
we shall now explore in greater detail. 



A. Coherent to incoherent transition 

We now return our attention to the crossover from 
coherent to incoherent transfer, defined by Eq. (|28[). 
Intuitively, we might expect the dynamics in the low- 
temperature (or weak-coupling) regime to be coherent; 
for example, in Fig. Q] incoherent relaxation only occurs 
in the high-temperature limit. If we therefore assume 
that the crossover itself occurs in the high-temperature 
regime, it is possible to derive an analytic expression gov- 
erning the crossover temperature by approximating the 
rates T y and T z . Details of this approximation, and its 
range of validity, can be found in Appendix [X] Gener- 
ally, for high enough temperatures and/or strong enough 
system-bath coupling (such that j3Vu <C 1) we can ap- 



proximate j xx (77) « J yy (r]) sa J yy (0) in T y and Y z , where 



7 TO (0) 



jj^ e 't>oCo(x,y) 



(30) 



with 4>o — 2tt 2 a / ljq f3 2 , x = nd/cf3 and y = u> c f5. The 



functions Co(x,y) and C2(x,y) are given by Eqs. (IA7|) 
and l|A8p , and the renormalisation factor B by the prod- 
uct of Eqs. (|A14I) and (|A15[) . If we further assume 
that the energy shift A3 vanishes in the high-temperature 
limit, Eq. (|28|) reduces to 



(r 2 - r„) = w R , 

and we arrive at the expression 



k B T 

UJQ 



V B e <l>oCa(x,y) 



w 4y/2ir 3 aC 2 (x, y) 



(31) 



(32) 



with solution, T c , giving the crossover temperature sep- 
arating the coherent and incoherent regimes. 

The dependence of T c on the various parameters in- 
volved in the problem is not straightforward, owing to 
the temperature dependence in the renormalisation fac- 
tor B, in the functions Co and C2, and in cf>Q. In 
fact, there are three distinct and important tempera- 
ture scales which determine when coherent or incoherent 
processes dominate: Tq = wo/(v / 2cv7r^s), which depends 
upon the system-bath coupling strength; T x = c/ dirks, 
which arises due to the fluctuation correlations and be- 
comes unimportant in the uncorrelated case (T x — > as 
d — > 00); and T y — u) c /k,B, dependent upon the cut-off 
frequency, and irrelevant in the scaling limit (y — > 00). 
Hence, changes in any of a, d, or lo c can have an effect 
on the crossover temperature. For example, the main 
part of Fig. [3] shows the solution to Eq. (|32t . i.e. the 
crosover temperature T c , as a function of the dimen- 
sionlcss cut-off frequency uj c /ujq. A calculation using 
Eq. (|28[) with the full rates, and including A3, is also 
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FIG. 3: Crossover temperature separating the coherent and 
incoherent regimes against cut-off frequency, for levels of cor- 
relation given by jj, — c/uiod = 0, (X = 0.5 and (x = 1, in- 
creasing as shown. The solid blue curves have been calcu- 
lated from Eq. (|28l) (using the full rates), while the dashed 
red curves are solutions to the high-temperature approxima- 
tion, Eq. (|32p . The inset shows the dependence on the level 
of correlation = oJod/c) for different cutoffs, uj c /ujo = 2, 
ui c /ujo — 3, and uj c /uiq — 4, again increasing as shown. Pa- 
rameters: a = 0.05 and V/uo = 0.5. 



shown for comparison. The three pairs of curves corre- 
spond to increasing levels of correlation, ordered as indi- 
cated. We see that, except for small lOc/^Jq m the case 
fx = 0, where A3 becomes important, solutions to Eq. (|32|) 
give an excellent approximation to the crossover temper- 
ature calculated using the full rates. This confirms that 
the coherent-incoherent crossover does indeed occur in 
the high-temperature (multiphonon) regime, and conse- 
quently could not be captured by a weak system-bath 
coupling treatment. 

As the cut-off frequency is increased from its minimum 
value, the crossover temperature begins to decrease. This 
behaviour can be understood qualitatively by examining 
Eq. pip , and considering the competition this condition 
captures between the rate T z — T y and the coherent inter- 
action Vr in denning the nature of the dynamics. Larger 
values of the cut-off frequency correspond to smaller val- 
ues of the renormalised interaction strength Vr (see e.g. 
Eq. (|A14p ). while the rates T y and T z vary less strongly 
with uj c in this regime. Thus, increasing lo c from its mini- 
mum value decreases Vr, and therefore reduces the range 
of temperatures for which 4Vr > T z — T y and coherent 
transfer can take place. Thus, the crossover temperature 
falls. Physically, this can be understood by noting that 
as the cut-off frequency is increased, so too is the effec- 
tive frequency range and peak magnitude of the system- 
bath interaction, characterised by the spectral density 
[Eq. Hence, increasing from small uj c /ujo, the en- 

vironment begins to exert an enhanced influence on the 
system behaviour, and so coherent dynamics no longer 
survives to such high temperatures. As uj c continues to 
increase, however, we see the crossover temperature then 
begins to rise. The renormalisation factor B tends to 



zero with increasing uj c and here becomes the dominat- 
ing quantity, thus causing the rate T z — T y ~ 0(B 2 ) to 
vanish faster than the renormalised donor-acceptor cou- 
pling V R = BV. 

The interplay between the size of lo c and the level of 
spatial correlation is best understood by considering the 
inset of Fig. [3] For all curves shown the crossover tem- 
perature increases as the distance d is reduced, since the 
level of correlation fi increases correspondingly. As we 
have seen previously in Fig. [21 stronger correlations al- 
low coherent dynamics to be observed at higher temper- 
atures; since environmental effects arc suppressed, so the 
crossover temperature T c must rise. This behaviour can 
be attributed to an increase in the renormalised interac- 
tion strength, Vr, in relation to the rate V Z —T yi this time 
with variations in the correlation level /i. Interestingly, 
as the cut-off frequency is increased up to w c /wo = 4 
(lowest curve), we see that not only does the crossover 
temperature decrease, but also that the degree of corre- 
lation necessary to show a marked rise in T c increases. 
As can be seen by comparing the separation between the 
different curves in the main part of the figure, increas- 
ing the cut-off frequency tends to suppress the extent to 
which correlations are able to protect coherence in the 
system. This tallies with the dynamics shown in Fig. [21 
for which u; c /cjo = 4, and correlations as high as /i = 2 
were needed before a significant change in behaviour was 
seen. Finally, since the renormalisation factor B tends 
to a constant non-zero value as the correlations vanish at 
large d (as opposed to B — > as uj c — > 00), the depen- 
dence of the crossover temperature on /1 is monotonic, in 
contrast to its dependence on uj c . 



IV. OFF RESONANCE 

It is often the case in practice that the donor and ac- 
ceptor will have different excited state energies, 61 — 62 = 
e^0, and so we now turn our attention to energy transfer 
dynamics under off-resonant conditions. Regarding the 
coherent to incoherent transition, in the resonant case we 
were able to identify this point with a pair of conjugate 
eigenvalues converging on the real axis, thus changing 
oscillatory terms into relaxation. We might hope that 
in the off-resonant case we are able to establish a simi- 
lar crossover criterion, and again use this to investigate 
the effects of bath correlations and the cut-off frequency. 
However, we shall see that such an identification is less 
straightforward in the off-resonant regime. 

We first present the full Bloch equations describing the 
evolution of our donor-acceptor pair for arbitrary energy 
mismatch. As in the resonant case, we have an equation 
of motion of the form a. = M ■ a + b, but now the matrix 
M is given by 



M 




-(e + Ai) 

-r„ -2BV R 
b- 1 (2v r + \ 3 ) -r. 



(33) 



8 




20 40 60 80 100 

U) \ 



Fig. Q] (plotted in the resonant case) and Fig. [4] we see 
that the presence of a substantial energy mismatch causes 
the low-temperature population oscillations to increase in 
frequency but decrease markedly in amplitude, such that 
for fcflT/wo = 5 oscillations are now almost impercepti- 
ble. We also see that the population difference tends to a 
non-zero steady-state at low temperatures, as we might 
expect from simple thermodynamic arguments, since the 
states a z = 1 and a z = — 1 now have different energies. 
As the temperature is raised, however, the dynamics still 
looks to be approaching that shown in the resonant case 
of Fig. [I] 



FIG. 4: Population difference for an off-resonant donor- 
acceptor pair as a function of scaled time coot. Temperatures 
kBT/uoo — 1 (blue dashed curve), fcsT/uo = 5 (green dotted 
curve), ksT/uia — 12 (orange solid curve) and fesT/wo = 20 
(red dot-dashed curve) are shown. Parameters: a = 0.05, 
V/wo = 0.5, e/wo = 1, Uc/^o = 4, and /i = 0.5. 



with b = (—Bk x , —Bk v , —k z ) t . The rates become 



(34) 



T, = 2V 



{ W 2 e 2 \ 

( -fr«(o) + (r/xx(v) + 7x*(-*?)) ) 



(35) 



with T z — T x + T y , and the energy shifts 

2V 2 e 

Al = —^-(Syy{ri) ~ Syy(~T])), (36) 

2V 2 e 

A 2 = (S xx (r]) - S xx (-r])), (37) 

V 

W 2 V Ri 



A. 



V 



(Syy(v) ~ Syy(-V)). (38) 



The remaining quantities are 

4V 2 Vre ( 1 \ 

C = ( 7**(0) ~ -fil*M + 7^(-r?))J , (39) 

2V 2 V R 

K x = (jyyiv) ~ lyy{-ri)), (40) 

8V 2 V R e ( 1 \ 

Ky = \s*M - ^{s x M + s xx (- v ))) , (4i) 

V 2 e 

K z = {{lxx{v) - lxx{~V)) + (lyyiv) - lyy(-V))) ■ 

(42) 



Here, i] — ^/ e 2 + is the system Hamiltonian eigen- 
state splitting in the polaron frame. 

To exemplify the dynamics generated by the full off- 
resonant Bloch equations, in Fig. 2] we plot the evolu- 
tion of the population difference in the case of donor- 
acceptor energy mismatch, e = 2V. By comparison of 



A. Correlated fluctuations 

Let us now look at the effect of correlated fluctuations 
in the off-resonant dynamics. Though analysis of the 
full Bloch equations is now more complicated than in 
the resonant limit, we should still expect changes in the 
level of donor-acceptor fluctuation correlation to have a 
qualitatively similar effect on the transfer process as out- 
lined in Section [TTTT To illustrate that this is indeed the 
case, in Fig. [5] we plot the donor-acceptor population dy- 
namics under off-resonant conditions for three different 
levels of fluctuation correlation (increasing from top to 
bottom) . Just as we found in the resonant case of Fig. [5J 
an increase in correlations enhances the lifetime of co- 
herence present in the energy transfer process, and can 
even move the dynamics from a high temperature (or 
strong-coupling) predominantly incoherent regime to an 
effective low temperature (or weak-coupling) regime dis- 
playing pronounced coherent oscillations. In addition, 
in the off-resonant case stronger correlations also serve 
to amplify the coherent contribution to the full energy 
transfer dynamics (made up of distinct coherent and in- 
coherent parts, as we shall show below), since the renor- 
malised interaction strength Vr increases in relation to 
the energy mismatch e. 

In an effort to put these qualitative observations on a 
more quantitive footing we could analyse the eigensystem 
of the full off-resonant M [Eq. (13"3"1) ]. in a similar manner 
to the resonant case. However, finding the eigensystem of 
M is now far less straightforward and analytical solutions 
to the full Bloch equations are consequently lengthy, and 
therefore of little direct use in gaining an understanding 
of the behaviour seen in Figs. |4]and[5] The rest of this 
section is thus devoted to deriving simplified expressions 
for the energy transfer dynamics in two important limits: 
(i) weak-coupling, or strong correlations, where coherent 
dynamics can dominate, and (ii) high temperatures (and 
weak correlations) , where the dynamics is similar in both 
the resonant and off-resonant cases. These expressions 
not only provide insight into the off-resonant behaviour 
of the system and the effect of correlated fluctuations, 
but also serve to highlight the difficultly in now defining a 
simple crossover criterion, as was possible in the resonant 
case. 
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FIG. 5: Population difference as a function of scaled time 
cjoi for temperatures of /cbT/cjo = 5 (blue dashed curve) and 
fcsT/wo = 10 (red dotted curve), and for separations corre- 
sponding to no fluctuation correlations, [i — c/uod — (top), 
weak correlations, /i = 0.5 (middle) and strong correlations 
fi = 2 (botttom). The insets show the evolution of the corre- 
sponding coherence a y . Parameters: a = 0.05, V/ujq = 0.5, 
e/ujQ — 0.5, and w c /uio = 4. 



B. Weak-coupling (or strong correlation) limit 

We begin by considering the weak system-bath cou- 
pling limit, which we obtain by expanding all relevant 
quantities to first order in J{uj). In fact, for strong 
enough fluctuation correlations this limit is attainable 
even if the system-bath coupling is not weak and/or the 
temperature not low, due to the factor (1 — Fr>(uj,d)) 
appearing in Eqs. (J7]) and (JUJ). With reference to our 
expressions for the correlation functions [Eqs. (|10p and 
pip ] , we see that within this approximation A xx (r) — > 



remains finite. We may then set to zero 
all rates and energy shifts which are functions of A xx {t) 
only in Eq. (1331) . This results in the far simpler form 




where the weak-coupling rate is given bv2i 



(43) 



T w = An\-^j J{v)(l~ fb(7,,d))coth(/V2), (44) 

and the two energy shifts may be written Ai = 
(e/?7)A and A 3 = (2V R / V )A, with A = 2V 2 (S yy (r 1 ) - 
S vv (—rf)). The inhomogeneous term becomes bw = 
{— Bk x , 0, — (c/2Vr)k x } t in the same limit, which leads 
to the weak-coupling steady state values of 



a x {oo) 
a z (oo) 



2BV R 



tanh(/3??/2), 



V 



■tanh(/3?7/2), 



(45) 
(46) 



and a y (oo) — 0. As in the resonant case (in which 
there was no weak-coupling approximation) , this steady- 
state has the same form as that expected from a stan- 
dard weak-coupling approach, though with the replace- 
ment V — > Vr, and the extra factor of B suppress- 
ing the coherence a x (oo). As the energy mismatch in- 
creases in relation to V, the weak-coupling steady state 
therefore becomes increasingly localised in the lower en- 
ergy state 1 2) = \GX). Interestingly, this contrasts with 
the qualitatively incorrect form (at low temperatures at 
least) given by the Non-Interacting Blip Approximation 
(NIBA), a™ BA (oo) = -tanh09e/2),22jZ2 which predicts 
complete localisation in the lower energy state at zero 
temperature, regardless of the size of e/V. We should 
thus expect the present theory to fair far better than 
the NIBA for low-temperatures (or weak-coupling) in 
the off-resonant case, e ^ 0. The rate Tyy given in 
Eq. dHJ) is also of the form expected from a weak-coupling 
treatment, though once more with the renormalisation 
V — > Vr. In fact, such a replacement is sometimes 
made by hand in weak-coupling theories to provide agree- 
ment with numerics over a larger range of parameters^ 
though it arises naturally in the polaron formalism here. 
We can therefore conclude that, in addition to allowing 
for the exploration of multiphonon effects , n i 12 i 51 i 55 the 
polaron master equation provides a rigorous way to ex- 
plore the (single-phonon) weak-coupling regime for spec- 
tral densities of the type in Eq. (f29|) £!■ 

As before, to find the time evolution of a we evaluate 
the eigensystem of Mw and use Eq. (fT5f . For the initial 
state a(0) = {0,0, 1} T we find population dynamics 



a z {t) = 



-T w t 



- 1 



tanh(/3r?/2) 



— ^e 2 cos — — - - — sin — — , 
V 2 V ^ 2 i £ w V 2 J J' 



(47) 
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where the weak coupling oscillation frequency is given by 



£w = ^(77 + A) -r^, 



(48) 



which we expect to be real to be consistent with our 
original expansion. The first term in Eq. (1471) . propor- 
tional to (e/rf) and present nowhere in the resonant case, 
describes incoherent relaxation towards the steady state 
value given by Eq. (l46|) . The second term, proportional 
to (Vfl/77) 2 and having a similar form to the resonant 
dynamics, describes damped oscillations with frequency 
£w Importantly, these oscillations have a temporal max- 
imum amplitude of iV^/r] 2 < 1, compared to 1 in the 
resonant case. The effect of the energy mismatch in this 
limit is thus to suppress the amplitude of any oscillations 
in the population difference, while increasing their fre- 
quency due to the dependence of on V in Eq- (l48lh 
exactly as seen in Fig. 



C. High temperature (or far from resonance) limit 

At high temperatures and weak correlations, we find 
that the population dynamics appears to be relatively in- 
sensitive to the size of the energy mismatch. In order to 
investigate this effect in more detail, we shall now make 
a high-temperature (or strong system-bath coupling) ap- 
proximation to the full energy transfer dynamics. 

Specifically, we consider the regime Vr/c <C 1. This 
limit can in fact be achieved in two possible ways. Firstly, 
recalling that Vr = BV, we see that Vr can be made 
small by increasing the system-bath coupling strength 
or temperature, such that B <C 1. Alternatively, if the 
donor-acceptor pair are far from resonance, the ratio V/e 
will be small, and hence Vr/e smaller still. Observing 
that the correlation functions given by Eqs. (ITU1) and (fTTj) 
are both proportional to B 2 , we can see that all dissipa- 
tive terms in the equation of motion, a = M ■ a + b, 
are at least of order Vr. We proceed by keeping only 
terms up to order (Vfj/e) 2 in the full off-resonant M and 
b. This allows us to set A3, £, n x and n y to zero, while 
the remaining quantities reduce to 



(49) 



T z = V 2 (j xx (ri) +'y xx {-r)) +73,3,(77) +l m {-rf)), (50) 
\ 1 =2V 2 (S yy (r ] )-S yy (-r 1 )), (51) 
A 2 = 2V 2 {S xx (i 1 )-S xx (- V )) 7 (52) 

Kz = V 2 (j xx (lj) - J xx (-T]) + 7j,j, (7?) - lyy{-V))- (53) 

Hence, in the high-temperature limit, Eq. (I33[) takes on 
the simpler form 




(54) 



while the inhomogeneous term reduces to bjjT = 
{0, 0, — k z } t . We then find the approximate steady-state 



population difference 

a z (oo) = 



1 + ^W^-l) )ianh( 0//2). (.V,) 



valid up to second order in Vr/c. For Vr <C e, this 
steady-state is strongly localised in the low energy state 
(a z (oo) w —1) if e ^ fc^T, though for e <C ksT thermal 
effects dominate and a z (oo) « as in the resonant case. 
Again, this behaviour tallies with Fig. @] 

To obtain the corresponding population dynamics, we 
note in reference to Eq. (fT5|) that the coefficients a^, the 
eigenvectors mi , and the eigenvalues qi will contain pow- 
ers of our expansion parameter Vr/6. Expanding both qi 
and the products aiirii to second order, we find 



a z {t) 



-r 2 t 



-(1 



1 



4Vg 



^c- r ^ 2 cos(ei) 



) tanh 



4^1 



r. 



- i 



(56) 



where the shifted oscillation frequency is 

e = e + (l/2)(Ai + A 2 ) + 2e{V R /e) 2 . (57) 

As in the weak-coupling case [Eq. (|47p] the evolution of 
the donor-acceptor population difference consists of two 
contributions; incoherent relaxation towards the steady- 
state, and an oscillatory component with vanishing am- 
plitude as Vr/c — > 0. The energy mismatch again serves 
to suppress oscillations in the population difference. 

The most striking feature, however, of Eq. (|56|) is that 
there is an oscillatory component at frequency e at all. 
In the high-temperature limit, we might expect that this 
frequency would reach a point where it becomes imagi- 
nary and a z (t) displays purely incoherent relaxation, as 
in the equivalent resonant case. However, we can see 
that this is not the case since e is always real by defini- 
tion. Furthermore, at very high temperatures e — > e, and 
it therefore also remains finite. Eq. (|56[) thus highlights 
an important difference between the energy transfer dy- 
namics in resonant and off-resonant situations. In the 
resonant case, as temperature is increased, the energy 
transfer process becomes less coherent through a reduc- 
tion in oscillation frequency (i.e. Vr becomes small in 
comparison to T z — r y ), eventually reaching a point at 
which population relaxes incoherently towards the steady 
state. In the off-resonant case, the transfer process be- 
comes less coherent predominately through a reduction in 
oscillation amplitude. For high temperatures, an oscilla- 
tory component is still (in theory) present in the system, 
although it becomes ever more dominated by incoherent 
relaxation towards the steady-state population distribu- 
tion, which depends upon the ratio e/ksT. These fea- 
tures are clearly seen in Fig. |4j 

Only to first order in Vr/c do our expressions predict 
purely incoherent off-resonant population transfer: 



a z {t) 



-r 2 t 



(1 



-T z t 



)tanh(/3r7/2). (58) 
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FIG. 6: Coherence {a y ) as a function of scaled time toot for 
resonant (e = 0, solid curve) and off- resonant (e/wo = 0.2, 
dashed curve) cases. The temperature, feT/wo = 13, is cho- 
sen to be above the relevant crossover T c in the resonant case, 
such that the resonant dynamics is guaranteed to be incoher- 
ent. Parameters: a — 0.05, V/ljo — 0.5, w c /wo = 4, and 
H — 0.5. The inset shows the corresponding population dy- 
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FIG. 7: Coherence (a y ) as a function of scaled time LO^t for 
temperatures of ksT/uio = 1 (blue dashed curve), ksT/ujo = 
5 (green dotted curve), fcsT/wo = 12 (orange solid curve) and 
UbT /uo — 20 (red dot-dashed curve). Parameters: a = 0.05, 
V/ljo = 0.5, cj c /uo — 4, e/cJo = 2 and /i = 0.5. The inset 
shows the corresponding population dynamics. 



Let us also consider the evolution of a y in the same limit: 
ay (t) = _^V(i/2)r,t sin( -^ (5g) 

Hence, although the donor-acceptor population itself 
evolves entirely incoherently in this limit, the coherences 
may still perform oscillations due to the energy mis- 
match. To illustrate the difference in the transition to 
incoherent population transfer on- and off-resonance, in 
the main part of Fig. [6] we plot the evolution of the co- 
herence a y (t) in both cases. The parameters have been 
chosen such that the resonant dynamics is in the incoher- 
ent regime (T > T c ), hence the resonant a y displays no 
oscillations [see Eq. (ffifll ]. In accordance with Eq. ([59]). 
however, the introduction of an energy mismatch induces 
oscillations in the donor-acceptor coherence. While these 
oscillations have an almost negligible amplitude, this be- 
haviour serves to illustrate the subtlety in defining a strict 
crossover from coherent to incoherent dynamics in the off- 
resonant case. In particular, despite the different forms 
of coherence behaviour, the corresponding (essentially 
incoherent) population dynamics shown in the inset is 
almost indistinguishable in the two cases, even though 
there should still be a strongly suppressed coherent con- 
tribution in the off-resonant curve. 

An alternative way to obtain oscillations of the coher- 
ence oty in a regime of predominantly incoherent popu- 
lation transfer is to introduce a large energy mismatch 
(i.e. make V/e small) at low temperature, as shown in 
Fig. [7] Here, for the lowest temperature considered the 
population relaxes towards its steady state value with 
little sign of oscillation, while the coherence performs os- 
cillations with a significant amplitude and considerable 
lifetime. This behaviour is strongly suppressed, however, 
as temperature increases, such that ksT > e. 



V. SUMMARY 

Motivated by recent experiments which suggest that 
quantum coherence can survive in energy transfer pro- 
cesses even under potentially adverse environmental con- 
ditions, 1-9 we have investigated various factors that de- 
termine the nature of the energy transfer dynamics in 
a model donor-acceptor pair. To do so, we used a po- 
laron transform, Markovian master equation technique. 11 
This formalism is attractive as it allows for exploration of 
both the low-temperature (or weak-coupling) and high- 
temperature (or strong-coupling) regimes, as well as reli- 
able interpolation between these two limits, provided the 
ratio V/oj c does not become too large i 13 ' 51 ' 55 We are also 
able to consistenly describe off-resonant effects, unlike in 
the NIBA ) 69 i 72 and the influence of bath correlations. 

In the resonant case we identified a crossover temper- 
ature separating coherent and incoherent energy trans- 
fer. We found a non-trivial dependence of this temper- 
ature on both the degree of spatial correlation within 
the bath-induced fluctuations, and also on the cut-off 
frequency of the bath spectral density. Smaller cut-off 
frequencies were found to enhance the extent to which 
bath spatial correlations are able to protect coherence in 
the system. The crossover generally occurs in a high- 
temperature limit where multiphonon effects dominate, 
and so could not be captured by a standard perturbative 
treatment of the system-bath interaction. 

In the off-resonant case we found that coherent and 
incoherent regimes are less easily defined. In particu- 
lar, for a sufficiently large energy mismatch between the 
donor and acceptor, coherence can in theory be present 
at all but infinite temperatures, albeit with an ever de- 
creasing amplitude. However, using analytic expressions 
derived in various limits, we were able to characterise the 
off-resonant energy transfer process over much of the pa- 
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rameter space, illustrating the suppression of coherence 
in the population dynamics with increasing temperature 
or energy mismatch. We also showed that strong corre- 
lations have a qualitatively similar effect to the resonant 
case, protecting coherence in the transfer process. 

While we have concentrated in this work on elucidat- 
ing general features of donor-acceptor energy transfer dy- 
namics using a simple model system, the insight we have 
gained could be relevant to a variety of systems. In addi- 
tion to those already mentioned^— closely-spaced pairs 
of semiconductor quantum dots could provide a solid- 
state implementation of the model studied here^ In par- 
ticular, our polaron master equation theory provides a 
bridge between the weak^L and strong^ system-bath cou- 
pling approximations already explored in this context. It 
would also be interesting to analyse the energy transfer 
dynamics of larger donor-acceptor complexes within the 
polaron formalism^ to see if further understanding of 
the interplay between coherent and incoherent processes 
in such systems could be obtained. Finally, it would 
be desirable to perform a thorough investigation of the 
regime of validity of the polaron approach by comparison 
to numerically exact techniques/® 
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Appendix A: High temperature rates 

Here we show how to obtain analytic approximations 
for the decoherence rates at high temperatures by use of 
a saddle point integration. Within our formalism there 
are two rates which need to be evaluated: 

n2 r+oo 

7xx {v) = =r drelT " (e0(r) + e " 0(r) - 2) ' (A1) 

* J-oo 
d2 r+oo 

(77) = — / dre lT,) (e 0(r) - e^ (r) ), (A2) 

* J -00 



Ivy 



where </>(r) is given by Eq. (fT2j) . With the appropriate 
manipulations 51 it is possible to show that 



7xx(r?) = ^e^/ 2 / dr<--"u-" lr '- <• -2i. 



(A3) 



B 2 



Ivviv) = ^r^' 2 / dre ir V W " e" 0(T) ), (A4) 
2 J -00 

where now <f>(r) = 4>{—t) = 4>{t — i/3/2), and is given 
explicitly in integral form by 

-. , f°° J[oS) , , cos(wr) , , , 

4>{t)=2 dtJ-V(l-F p (a;,d)) . ^ ' (A5) 







sinh(^w/2)' 



Using a super-Ohmic form of spectral density, J(w) = 
ctL0 3 L0Q 2 e~ u '/ u ' c , and assuming system-bath coupling in 
three dimensions such that Fd(lj, d) = sinc(wd/c), al- 
lows 4>{t) to be found analytically. We find </)(t) = 
4> C(x,y,T'), where 

C(x,y, T ') = 

1 1 i . , 
- + r'+: 

2 V 7T 



—I 

2lTX 



4> 



1> 



-(r'-x) 



1 1 1 i , , , 

+ i>\n + - + ~( T ~ x ) 

2 V 7T 



1 1 i . 

- + - + -/ + x) 

I y 7r 



1 



Here, <f>o = 2n 2 a/(uj 2 (3 2 ), x — nd/c(3, y — uj c (3, t' = 
7tt//3, tp(z) is the the digamma function, and ip'(z) its 
first derivative. 

To proceed, we assume a high-temperature or strong- 
coupling regime, such that the dominant contributions 
to the integrals in Eqs. (IA3|) and (|A4I) will come from 
the peak in </>(r) at r = 0. More specifically, inspec- 
tion of C(x,y,i J ) reveals that for y » 1 (the scaling 
limit of large w c ), we require (f>o 1 for large x (weak 
correlations), or 4>qx 2 ^> 1 for small x (strong correla- 
tions), in order for an expansion of </>(r) around r = to 
be valid. These definitions of the high-temperature (or 
strong-coupling) regime tally with the expansion param- 
eters identified in Ref. [ll|. In the opposite limit, y <C 1, 
we generally need 4>ox 2 y 3 /ir 4 ^> 1, except in the limit of 
very large separations (vanishing correlations), x — > 00, 
where 4>oy 3> 1 is the relevant condition. 

With these conditions in mind, we therefore expand 
4>(t) to second order in r' = ttt//3, which gives 



(r)^MC (x,y)-r' 2 C 2 (x,y)), 



(A6) 



where 



Co{x,y) 



+ 



1 1 ix 
- + - + - 

2 y tt 

-%i> (0 + -)' 

TT Z V 2 V / 



1 1 ix 

2 y tt 



(A7) 



and 



C2(x,y) 



2tt 3 x 



1 1 ix 
- + - + - 

2 y tt 

1 



1 1 ix 

2 y tt 



(A8) 

In this limit, Eqs. (|A3|) and (|A4|) are dominated by the 
terms containing a factor of exp[</>(r)], allowing us to 
write 

0+00 



2tt 



(A9) 
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The integral is Gaussian and we arrive at the result 

(A10) 



2\/ 7r6 2(2:,2/)</>o 



which reduces to 



2 v / irC 2 (x,y)(t)o' 



(All) 



if the temperature is high enough such that l/rj/3 3> 1. 

It remains now to determine the bath renormalisation 
factor B, given by Eq. 0. To do so it is helpful to 
separate vacuum and thermal contributions. We write 
B = B B t h, where 



Bq = exp 

and 

B t h = exp 



J(w) 



(l-F D (u,d)) 



(A12) 



do; 



J(w) 



(l-F D (w,d))(coth( i Sw/2)-l) 
(A13) 



Inserting the spectral density, and again assuming 
system-bath coupling in three dimensions, we find 



Bo = exp 

and 

Bth = exp 



Wq \1 + (du) c /cy 



(A14) 



ix/n)—H(y 1 +ix/it)) 



2^'(1 + V- 1 ) 



(A15) 



where if (to) — ^ —1 (1 /i) is the mth harmonic number. 
We note that in the infinite separation (uncorrelated) 
limit, one finds B${d — > oo) = exp[— a^u^/wo)], and 

Bth = exp[a( Wc 2 /^ 2 )(l - jT 2 (C(2, 1 + y' 1 ) + C(2, iT 1 )))], 

(A16) 

where C( s i a ) is the generalised Riemann zeta function. 
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